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ABSTRACT 



When a differentially rotating magnetic star undergoes "turn-over", i.e., its 
magnetic symmetry axis is inclining at a gradually increasing angle with respect 
to its invariant angular momentum axis up to the perpendicular position (sit- 
uation known as "perpendicular rotator"), the poloidal magnetic field causes 
interchange of angular momentum via hydromagnetic Alfven waves propagating 
along the poloidal magnetic field lines and, accordingly, the resulting perpendic- 
ular rotator obtains uniform rotation. This dissipative mechanism is represented 
by the so-called "turn-over viscosity" . For several white dwarf and neutron star 
models, we compute the energy dissipation due to the turn-over viscosity; we 
also compare this effect with the corresponding dissipative effects induced by the 
electron, ion, and turbulent viscosity. 



Subject headings: methods: numerical — stars: magnetic fields — stars: neutron 
— stars: rotation — white dwarfs 



Introduction 



According to the "turn-over scenario" (for details see |Geroyannis 2001] , hereafter G01 



§ 3), the rotating magnetic star undergoes an axisymmetric "early evolutionary phase", 
during which the combined effect of rotation and (toroidal plus poloidal) magnetic field 
gives configurations with moment of inertia along the principal axis ^3 (coinciding with 
the rotation axis, the magnetic symmetry axis, and the invariant angular momentum axis, 
taken to be the z axis), J33, greater than the moments of inertia along the principal axes 
l\ and I2, In = I22', i-e., J33 > In. At some particular time, however, which signals the 
beginning of a "late evolutionary phase" , secular angular momentum loss due to the surface 
poloidal magnetic field, B s , induces configurations with I 33 (L X ) = I n (L x ) (relation valid 
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for a particular angular momentum L x ) and, after that particular time, configurations with 
I^{L) < hi{L) when L < L x . Thus, a "dynamical asymmetry" is established in the sense 
that the moment of inertia along the rotation axis becomes less than the moments of inertia 
along the other two principal axes. This inequality is due to the action of the toroidal 
magnetic field, which tends to derive prolate configurations; while rotation and poloidal 
magnetic field tend to derive oblate configurations. 

Dynamically asymmetric configurations tend to turn over spontaneously and thus to 
rotate along axis with moment of inertia greater than J33. So, they become "oblique rota- 
tors" ; remark that "perpendicular rotator" is an oblique rotator with "turn-over angle" (or 
"obliquity angle": the angle of the magnetic symmetry axis with respect to the invariant 
angular momentum axis), x, equal to 90° (the internal dynamics of an oblique rotator has 
been studied in detail by |Mestel fc Takhar 1972j , hereafter MT72). 

The "terminal" nonaxisymmetric rigidly rotating model becomes perpendicular rotator 
on a "turn-over timescale", txov, comparable with the "period timescale", tp = P n0 w/Pnow', 
i.e., t TO v ~ tp. Although any intermediate oblique model with \ < 90° undergoes a com- 
plicated rotation along axis passing through the center of mass but not coinciding with one 
of its principal axes (MT72, §§ 2-3), the terminal model rotates along its I\ principal axis, 
which axis is now coinciding with the invariant angular momentum axis taken, that is, the z 
axis. Apparently, its angular momentum vector remains equal to that of the starting model 
(MT72, §§ 1-2), while its angular velocity vector obtains same direction with that of the 
starting model (coinciding with the angular momentum direction) but different magnitude. 

During the turn-over phase, the turn-over angle increases spontaneously up to 90° on 
timescale txovi since the rotational kinetic energy of the model decreases from a higher level 
with turn-over angle \ ^ 0° to a lower level with turn-over angle x — 90°. At this level, the 
model reaches the state of least energy consistent with its prescribed angular momentum 
and magnetic field. Accordingly, the excess energy due to the differential rotation field of 
motions, defined by the angular velocity component S7 3 along the (spontaneously turning 
over) J 3 principal axis (coinciding, in turn, with the magnetic symmetry axis), is dissipated 
down to zero under the action of turbulent viscosity in the convective regions of the model 
(MT72, §§ 3-4). Thus, differential rotation along the -Z3 principal axis decreases down to 
zero and all of its effects disappear. 

On the other hand, it is difficult for an oblique model to sustain differential rotation 
along its I\ principal axis mainly due to the destructive action of the poloidal magnetic 
field (see, e.g., |Moss, Mestel, fc Tayler T990| , hereafter MMT90, §§ 1-3; also, |Moss 1992 



hereafter M92, §§ 1, 5-6). Indeed, when a model has magnetic field symmetry axis inclined 
with respect to the rotation axis, nonuniform rotation in the electrically conducting stellar 
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material shears the field and generates magnetic torques. Thus, there is a competition 
between the efforts of the magnetic stresses to remove rotational nonuniformities, and those 
of the rotational velocities to bury and destroy magnetic flux. If the magnetic field and the 
electrical conductivity have appropriate values (see MMT90, § 4; also M92, § 6), then the 
magnetic field prevails and removes the nonuniformities of rotation. Hence, the terminal 
model rotates rigidly along its I\ principal axis with angular velocity flj.. 



2. The turn-over viscosity 

The turn-over timescale, tpov-, is given by (G01, eq. [41]) 



VEdw , n .- a / ff 2 \ '/l t V,\-'W\' ... 
trov^ — m (— 2 ) (__)^_\jA t )„ , (1) 

where the subscript sz denotes that averages calculated over the surface zone with base at 
the transition layer, £ s , and top at the boundary, £ 1; of the star (G01, eqs. [13]-[18]); averages 
without subscripts denote global averages. This trov estimate results as the mean energy 
per unit volume available for dissipation, DEdrd/V (DEdrd is the excess rotational kinetic 
energy due to differential rotation and V is the volume of the configuration), divided by the 
mean energy dissipation during the turn-over due to turbulent viscosity in the convective 
surface zone per unit volume and per unit time (i.e., the product of the last 5 terms in the 
right-hand side of eq. for details see MT72, eqs. [7], [49]-[50]). The most quantities in 
equation (JT|) are calculated under the assumption of slightly declining from sphericity. In 
particular, the density is approximated by g(£, v) ~ £?(£)> the pressure is approximated by 
P(£, v) ~ -P(O) the local sound speed, a, is approximated by 



P(() (2) 



M £>«) ' 

the turbulent mixing length, A t , results from the "mixing- length theory" (MLT; see, e.g., 
Uanuto fc Mazzitelli 1991 , hereafter CM91) as A t = A P h P (CM91, §§ 2.8, 3.2), where the 



fine-tuning parameter Ap is taken equal to unity and hp is the pressure scaleheight, 

A t (t v) ~ A t {£) ~ h P (0 = A MLT -52- = A MLT ; (3) 

A-mlt is the standard MLT fine-tuning parameter taken to be of order 10 _1 , M(£) is the 
mass inside a sphere of radius £, and the turbulent velocity, V t , is approximated by (see, e.g., 
Brickhill 1990| , hereafter B90, eq. [5]; in combination with CM91, eq. [38]) 



\ 0(0 ■ { ) 
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An interesting question is what kind of globally acting "turn-over viscosity", fiTOV, 
can substitute the locally acting "turbulent viscosity", fi t , so as to induce energy dissipated 
over the time t TO v, DE TOV) equal to the difference DE DRD = T xx — T RR of the rotational 
kinetic energies of the starting and the terminal models. The reason for seeking a global 
viscosity, instead of the local turbulent viscosity, is that differential rotation, which causes 
this kind of viscosity, is a global property of the model. We can first substitute the local 
turbulent velocity field, V t , by a global "turn-over velocity field", Vtov ~ R/trov ~ Knp 
(the latter symbol denotes the mean speed of the magnetic poles on trajectories coinciding 
with circumference quadrants). Then the turbulent viscosity (see, e.g., B90, eq. [8]) 

MO = t^MO^MO ~ (V t ) sz (h P ) sz (5) 

can be substituted by the turn-over viscosity 

/WO * ~<?(0 (Vtov) A h TOV (h P ) ~ ^(O^g^, (6) 

(since for our models (hp) ~ R/10 2 ) where A TOV is a fine-tuning parameter. Furthermore, 
the energy dissipated due to such turn-over viscous friction per unit time, Dtov, is given by 
(see, e.g., |Geroyannis fc Sidiras 1995| , eq. [7]) 

D TOV = 4vra 3 ^ c J 1 ^ TOV ^) ^ (1 _ ^ t*A ^ (?) 

and assuming that Dtov remains constant over the time trov the energy, DEtov, dissipated 
due to such viscous friction over the time tTov is equal to t T ovDTov', so, this equality leads 
to the determination of the fine-tuning parameter A TOV . 

A second step then is to develop a phenomenological model involving further details 
on the time evolution of the turn-over process (G01, § 6). In the so-called "mixer-mixture 
model" (MMM), we assume that the turning over configuration acts as a "mixer" on the 
"mixture" configuration, which remains almost axisymmetric with respect to the invariant 
angular momentum axis, i.e., the z axis (the meaning of "almost axisymmetric configuration" 
is that the angle 7 between its angular velocity and its invariant angular momentum remains 
small: 7 < 1°, say; apparently, 7 becomes zero after the termination of the turn-over; in 
the following, we shall assume that such a deviation of the mixture configuration from 
axisymmetry is negligible and we shall treat it as being axisymmetric with respect to the z 
axis). The mixer consumes its excess rotational kinetic energy due to differential rotation 
in the process of homogenizing the rotation of the mixture. In particular, when viewed in 
the mixture's system of reference, the mixer generates a velocity field, which carries mass 
elements across cylinders of different angular velocities and momenta. This is due to the 
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fact that the rotational circular orbit of a mass element of the mixer, as seen in the mixer's 
system of reference, crosses several cylinders with different angular velocities, as seen in the 
mixture's system of reference. Thus, there is an interchange of angular momentum, which 
leads to uniform rotation of the mixture. 

For (1) the turn-over angle \ (with u x = cos(x)), (2) the central angular velocity fl^ of 
the mixer, and (3) the differential rotation reduction factor F r of the mixer, we adopt the 
time dependent scheme (G01, eqs. [52]-[54]) 

/ i + e (*-i)« \ /i + e i>-i)A 

QiQinut) = l 1+e ( K - l/L (t) )S J QoiQini) ~ ( j + ^ J Qo{Qini), (8) 

where the so-called "threshold parameter", k, is taken equal to k = cos(45°) ~ 0.707, the 
so-called "delta parameter", 5, is chosen equal to 5 ~ 5, and the auxiliary function ^(t) 
varies linearly with the time in the interval [0, troy] 5 

When 5 is relatively large (8 = 50, say), the second term in the right-hand side of equation 
(|) becomes negligible with respect to its counterpart, and QoiQini) can be simply set equal 
to the initial value(s) Q inh v x[ini] ^ 1, Q 3c [ini] = 0(L ra ), and F r[ini] = 1.00, or 0.50. However, 
as our numerical results show, the present study requires 8 ~ 5; since then the second term 
in the right-hand side of equation (||) cannot be neglected, Qo{Qini) must be normalized to 
the value 

QoiQini) — ~ / 1+e ( K -l)4\ > 

(10) 



n 



so that Q(Qini,0) = Qi 



In the framework of MMM, we use a simplifying description according to which the 
turn-over speed at a particular point is proportional to the relative difference of the linear 
speeds due to differential rotation, v mxr {£,v mxT ,t) and v mxt (^,iy mxt ,t), in the mixer and 
in the mixture — where the rotation of the latter results as superposition of a gradually 
established rigid rotation about the z axis with angular velocity Q z (which is not involved in 
the calculation of the turn-over speed) and of a differential rotation about the z axis identical 
to the differential rotation of the mixer about the ^3 principal axis; for simplicity, we shall 
use the same symbols for the central angular velocity fl^ c and the angular velocity field Q3 
for the differential rotations of both the mixer and the mixture, keeping in mind that the 
only difference concerns their rotation axes — (G01, eq. [55]), 
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Vtov = V? 



mp 



l^mxtyC ^mxt-i t) V mxr (£,V mxr ,t)\\ 

|w(s(f , Vmxt), F r (t))s(£, Umxt) ~ , iw), F r (t))s(£, V mxr )\ 



uj{s{£,v mxt ),F r (t))s{£,v mxt ) J 

(11) 

where V mp is the expected order of magnitude for Vrov- The relation between the coordinates 
v mxt and v mxr of a particular point has as follows 



, +V X Vmxr ~ sfl - Vlyjl - V 2 mxr ifX + ^mxr<90° 

Vmxt{Vmxr,V X ) ~ \ fZ J f ~ 2 „ U^J 

-^ X *W + y 1 - ^ V V mxr « X + # m:r r > 90 

Assuming for simplicity that the coordinate z/ (without subscript) refers to the mixer con- 
figuration, we can write for the turn-over viscosity (G01,eq. [58]) 

fiTovtt, v, t) ~ z\q(0^ovVtov(C, v, t)A h TOV h P {i) . (13) 



Then the energy dissipated over the time trov, DEtov, is given by (G01, eq. [59]) 

ftTOV 

= / D TOV {t) dt = Ana 6 x 
Jo 



rtrov 

DE TO v 



s=s{£,v) 



(14) 



DTOv{t) is now a function of time, while in the foregoing interpretation it was taken approx- 
imately constant over troy- The required equality DE T ov — DE DRD is again the way for 
determining the fine-tuning parameter A TOV . 

On the other hand the energy, DE t , dissipated over the time trov due to turbulent 
viscosity acting over a "turbulent zone" [£ t , £i] (lying, in turn, within the surface zone [£ s , £i]) 
is given by 

rtrov 



ftTOV 

DE t = / D t (t) dt = 4ira° x 



(15) 
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where D t {t) is the power dissipated due to turbulent viscosity. Consequently, the required 
equality DE TO y = DE t can be used for accurately determining the base £ 4 of the turbulent 
zone. In the following section, we shall use this method for several white dwarf and neutron 
star models. Actually, we shall determine the so-called "representative turbulent zone" (RTZ, 
with & = C,rtz), which corresponds to an axisymmetric model with constant differential 
rotation, F r = F r ^ n q, and angular velocity, f2 3c = f^cpm], suffering from power loss due to 
turbulent viscosity, D t , equal to the average value (over the time £tov) of the power loss due 
to turn-over, (Dtov)^ t na t; is, 

WfiL M T f M£,MK 4 (1-- 2 ) ( d ^ S f^ ) 2 d^dv — (D TOV ) t . (16) 

In order for this method to be reliable, however, the power dissipated due to electron and 
ion viscosity of degenerate matter ( Purisen 1973| , § 3, eqs. [4], [6], respectively), Deiv 
(calculated by eq. with Utov substituted by the electron-ion viscosity /i e +i), should be 
small in comparison with the former, D t 3> Deiv'-, otherwise, the latter cannot be screened 
by the former and thus becomes significant in the framework of TOV scenario. 



3. Some results 

In this paper, we present some preliminary results. In particular, Table 1 gives the most 
interesting quantities involved in the subject discussed above. Our computations concern 4 
models of compact stars. Two of them simulate white dwarfs; the corresponding computa- 
tions have been made by the so-called "complex plane iterative technique" (CIT; jOeroyanms 



fc Papasotiriou 2000| , and references therein). The other two models represent neutron stars; 
the corresponding computations should be considered as rough computations. An accurate 
numerical treatment of neutron stars is now in preparation. 

As a summary, it seems that the representative turbulent zone is ~ 1% of the cor- 
responding surface zone of these models; this fact shows the high efficiency of turbulent 
viscosity in dissipating the rotational kinetic energy due to differential rotation and, thus, 
leading the configuration from an unstable high energy level (aligned rotator) to a stable low 
energy level (perpendicular rotator) within a relatively small turn-over timescale. 



Acknowledgments. The research reported here was supported by the Research Com- 
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Table 1. Summary of calculations 



Parameter 



Delta parameter (eqs. [8]-[10]), <5 (dimensionless) 

Initial turn-over angle, Xini (arcdegrees) 

Angular momentum, L xx 

Average surface poloidal field, B s 

Magnetic flux, / 

TOV timescale, t-^ov (yr) 

Present TOV time, t now (yr) 

Present central period, P n ow 

Average spin-down time rate due to turn-over, 

Present turn-over angle, Xnow (arcdegrees) 

Average power loss due to turn-over, (Dxov)t 

Power loss due to nonturbulent viscosity, D^iv 

Percent SZ^ thickness, (l ~ is X 100 

Percent RTZ h thickness relative to SZ, il^iSfZ. x 100 

Si— 5s 



Modcl a & Mass b 



WD d , 0.60 


WD, 0.89 


NS° , 1.61 


NS, 2.39 


1.00(+00) f 


4.41(400) 


1.00(4-00) 


1.00(400) 


1.41(400) 


1.41(400) 


1.32(400) 


1.41(400) 


7.00(448) 


2.17(449) 


2.31(4-48) 


2.00(448) 


5.00(4-06) 


5.00(406) 


3.00(4-13) 


2.00(413) 


4.01(4-24) 


2.28(424) 


4.44(4-25) 


1.90(425) 


2.05(407) 


3.00(406) 


1.60(4-01) 


9.81(402) 


1.36(407) 


1.56(406) 


8.52(4-00) 


8.55(402) 


1.42(402) 


3.31(401) 


2.14(-03) 


2.55(-03) 


8.90(-14) 


1.68(-13) 


1.67(-12) 


3.07(-14) 


7.11(401) 


7.21(401) 


6.30(401) 


8.31(401) 


1.71(431) 


1.45(433) 


4.06(4-41) 


6.19(439) 


8.42(428) 


3.65(430) 


6.50(+38) 


4.16(438) 


4.76(400) 


2.24(400) 


4.27(-01) 


2.30(-01) 


5.32(-01) 


2.15(400) 


5.47(4-00) 


2.52(400) 



Parameters of the axisymmetric differentially rotating starting model (aligned rotator) 

Central period, P xx 1.04(402) 2.48(401) 1.69(-03) 1.72(-03) 

Rotational kinetic energy, T xx 1.47(447) 1.81(448) 3.06(4-51) 2.55(451) 

Moment of inertia along the h axis, I DR ,u(h, L xx ) 1.79(450) 1.39(450) 9.15(4-44) 8.45(444) 

Moment of inertia along the I 3 axis, I DR ,3 3 (h, L xx ) 1.74(450) 1.38(450) 9.05(4-44) 8.23(444) 

Average surface poloidal field, B s 5.00(406) 5.00(406) 3.00(+13) 2.00(413) 



Paprameters of the nonaxisymmetric rigidly rotating terminal model (perpendicular rotator) 



Period, P RR 1.62(402) 4.07(401) 2.54(-03) 2.67(-03) 

Rotational kinetic energy, T RR 1.36(447) 1.67(448) 2.86(4-51) 2.35(451) 

Moment of inertia along the h axis, I RRi n (h, L xx ) 1.80(450) 1.41(450) 9.31(4-44) 8.49(444) 

Moment of inertia along the J 3 axis, /33(h) 1.72(450) 1.28(450) 8.63(4-44) 8.11(444) 

Average surface poloidal field, B s 5.05(406) 5.27(406) 3.10(+13) 2.02(413) 



a Basic model: F r = 1.00; h = 0.08 

b In solar masses, Mq 

c In cgs units, unless stated otherwise 

d WD: white dwarf (Chandrasekhar's model), accurate computations 

e NS: neutron star (Ideal Neutron Gas equation of state), rough computations 

'Parenthesized numbers denote powers of 10 

g SZ: surface zone 

h RTZ: representative turbulent zone 
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